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Abstract. The control protocols of two types of finite dimensional quantum 
systems are proposed. The feasibility of each protocol is possible and an arbitrary 
target state can be achieved from initial state by a constant field. The control 
parameters which are time periods of interaction between systems and control 
fields in each cycles are connected with the probability amplitudes of target states 
via trigonometric functions and can be determined analytically. 
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1. Introduction 

Quantum control proposed by Huang et. al. in 1983 pQ is to drive a quantum system 
from an initial state to an arbitrary target state through its interaction with classical 
control fields or with a quantum accessor. It has attracted much attention of chemists, 
physicists and control scientists due to its application in physics and chemistry, 
especially in quantum information processing. Various notations in classical control 
theory were generalized to the quantum control, such as open and closed control, 
optimal control [2 , controllability [3J 0] |B] , feedback control [7] and so on. Coherent 
and incoherent (indirect) control schemes are proposed. In later case the system is 
controlled by its interaction with a quantum accessor which is controlled by classical 
fields El [lOl [11] . Typically, in the approach of quantum control, one should first 
model the controlled system and examine its controllability which is determined by 
the system Hamiltonian and interaction Hamiltonian with classical fields, and then 
design classical fields to stream the system to the given target state, which is referred 
to as the control protocol and is the issue we would like to address in this paper. Some 
works were proposed along this line, for example, using the Cartan decomposition of 
Lie groups [12] . 

In our previous paper [13] . we proposed an explicit control protocol of finite 
dimensional quantum system using time-dependent cosine classical field. In this paper, 
we will use the square pules to control the finite quantum systems. Advantages of 
using square pules is that the interaction Hamiltonian in each control cycle is time- 
independent and thus can be treated easily. Two types of finite systems, the one with 
equal energy gaps except the first one and another one with all distinct energy gaps 
are considered. The relationship between probability amplitudes of target states and 
control parameters, the width of pules, is analytically established. 
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This paper is organized as follows. In Sec. II, we formulate the controlled system 
and control scheme and investigated the controllability. In Sec. Ill, we present the 
control protocol of system with all distinct energy gaps and in Sec. IV we consider 
the system with equal energy gaps except one. We conclude in Sec.V. 

2. Control systems and controllability 

2.1. Control Systems and protocol 

For an iV-dimensional non-degenerate quantum system with eigen energy E n and 
corresponding eigenstates |n), we suppose that its Hamiltonian described by 

JY 

H = Y,E n \n)(n\- (1) 

71=1 

Without losing generality, we assume H is traceless, i.e. trH — 0. In this paper, 
just two different types of systems arouse our interest, the first one of which having 
all equal energy gaps except the first one, namely 

Ml M2 = M3 = • • • = (J>N-1> ( 2 ) 

where /ij = Ei + i — Ei is the energy gap, and the another one with all district energy 
gaps 

IM^H, i?j = l,2,---,N-l, (3) 

For later convenience, we call them the system I and system II, respectively. We define 
energy gaps 

hu>i = E i+1 -E u i = l,2,---,N 
for system I, and 

Huii = E i+1 - E h i = l,2,---,N 

for system II. 

The purpose of this paper is to develop a control scheme to drive the systems to 
an arbitrary target state from an initial state, using some independent classical fields 
fm{t). The total Hamiltonian of the system and control fields can be generally written 
as 

M 

H = H + Hj, H I= J2 fm(t)H m , (6) 

m=l 

where M is the number of independent classical fields. 

For the iV-dimensional systems considered in this paper, the total control process 
includes N — 1 cycles. In the m-th cycle, we first apply a classical field f m {t) = d m 
to control the system for time period r m , and then turn off the control field such that 
the system evaluates freely for a time period T' m , as showed in Fig[TJ 

For system I, It is easy to find that nuj m ^ huj n for m =/= n. Thus we let the 
coupling between the field and system only causing the transition between E m and 
Ei , which described by the interaction operator 

fl m =|m+l)(l| + |l)<m + l|. (7) 



- 1, 



(4) 
(5) 
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Figure 1. Control fields. 



For system II, as all 7wl; m are different, we only make the transition between 
level E m and E m+ i occurs in each cycle. Then the process described by interaction 
operator 



H m = \m + l)(m\ + \m)(m + l\ 



(8) 



So the control process includes N — 1 cycles. 

For both systems, As the previous discussion, two processes implement in each 
cycle: Step (1) apply the control field to control the system for a time period r m ; Step 
(2) turn off the control field and allow the system to evaluate for a time period r' m . 
We will see that the first step provides the real probability amplitude of the target 
state and the second step provides the relative phases actually. Therefore, the control 
field can be rewritten as 

^m; ^m— 1 — t — tjn—l I'm 

0, otherwise, 



fm{t) 



(9) 



where t m — J2T=i( T k + T k)- Those N — 1 control fields are independent in the sense 
that each f m (t) ^ in different time period. 

The whole control process can be equivalently regarded as control by one control 
field f(t) = X)m=i fm(t), where /(*) is shown in Fig.l. 



2.2. Complete controllability 

After roughly presentation of the control protocol, we need to examine the 
controllability of this control scheme, namely, to examine whether the Lie algebra 
generated by the skew-Hermitian operators iHo and iH m 

C = Gen{iff , iH m \m = 1, 2, N-l} (10) 
is su(7V) [3]. Here H m for system I and II are 

iJ m = |l)(m + l| + |m + l)(l|, (11) 

H m = \m)(m + 1| + \m + l)(m| (12) 
respectively. Or equivalently, iH ,iH m generate the Chevalley basis of su(^V) [H[S] 

ix n = i(\n)(n + l\ + \n+ l)(n|), 

iy n = \n)(n + l\-\n+ l)(n| , 

ih n = i(\n)(n\-\n + l}(n + l\), (13) 
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where n = 1,2, ...,7V — 1. In fact, it is enough to prove ix n e C (or iy n e £), as 
= ^[iH^iXn] and ih n = -[ix n ,iy n ]/2. 
For system I, it is obvious that 

*a;i=iffi=*(|l)<2| + |2><l|)e£. (14) 

Then 

iy 2 = [iH2,m} = |2)(3| - |3)(2| e C, 

ix 2 = —[[iH^iHiliHo] = i(\2){3\ + |3}(2|) e C. (15) 

Recursively, we have 

iy m = [[iH m ,iH m -i],iHo] ■ (16) 

So the system I is completely controllable. 
For System II, iH m itself 

iH m = i(\m)(m + l\ + \m + l)(m\), (17) 

is nothing but the generator ix m . Therefore, the system II is completely controllable. 

3. Control protocol of system I 

In this section we illustrate the detailed control protocol of System I. Suppose the 
system is initially on the ground state IV'o) = |1}- The system is driven to an arbitrary 
target state after N — 1 cycles. Relationship between the control parameters {r m , T' m } 
and probability amplitude of target states is explicitly established. 

From the protocol we provide, there are two processes implemented in each cycle: 
l.fleid on ; 2. field off. In the m — th cycle, We define 

fiw m = E m+ i — Ei (18) 

and let i(£ m+ i + E\) = i.e. reset the energy base point. Then in each cycle, the 
system can be regard as just having two energy levels. Using these two conditions, the 
free Hamilton takes the form: 

N 

H = J2 E n\n)(n\ 

n 

N 

= £i|l)(l| + £ m+ i|m+l)(m+l|+ E n\n)(n\ 

71=2 

= -^|l)(l| + ^|m + l)(m+l|+ J2 E ^)(n\ 

= ^M™)+ ]T E n \n)(n\ (19) 

n = 2 

and the interaction Hamilton is 

HW = d m (\m + 1)(1| + \l)(m + 1|) = d m o-i m l (20) 
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So the total Hamilton becomes 
H = H + H {m) 

hus. 



N 



J- E n \n)(n\+d m a% 



m) 



n = 2 



N 



(21) 



n = 2 



where H { c m) = ^ai m) + d m ai m) . 

As does not depend on t, the time evolution operator can be written as 



U {m} (t) 



-iHt/h 



exp 



it 



t 



N 



x exp 



HW+ E ^»l»>(' 

11=2 

in(Q m t/h)H^ 
E ^»|n)<n|t/R 



cos(O m t/7i) — ^— sin 



71=2 



from the fact 



2n 



n 2n T 

7TL ' 



2n+l 



O 2 ™ FT C m ) 
iL m 1J -c i 



(22) 



(23) 



where fi m = ^(^) 2 + ^. 

For convenience, it won't change the hnal results naturally, if we let h = 1. So 



C/ (m) (<) = 



cos(Q m t) - — sin(a„*)#i m) 

1 Im 



x exp 



N 



-i E E n \n)(n\t 



n = 2 



(m > 2) 



(24) 



further more, the evolution of system in each cycle can be investigated in schrodinger 
picture, which is more convenient than the method using a changing field. 



3.1. Cycle 1 

In the first cycle, the initial state of the system is |1). Noticing that ([24| exclude the 
case m = 1, we has to calculate E/W(t) in cycle- 1, written as 

N 



cos(fiif) - — sin(fiif)i? P 
Hi 



exp 



\^E n \n){n\t 



n— 3 



.(25) 



Then after time period t\ in step one, system's state becomes 

hM =f/ (1) (n)|i) 



co S (Om) |1) + Bin(n m ) |1) - ^sm(fiiri) |2> (26) 
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when di 3> i-e. the magnitude of external field is very large. Then Tiuii/Qi — > 
0, di/Qi — > 1, we get 

=cos(n 1 r 1 )|l)-isin(n 1 r 1 )|2). (27) 
Therefore after free evolution for r{ time period, the final state is 

\^=e- iH ^\^)=a^\l)+a^\2), (28) 



with 



a[ 1} = cos(fi 1 ri) e - lEl ^ (29) 
4 X) = -isin^irOe-'^ 2 ^. (30) 



3.2. Cycle 2 



For cycle-2, the initial state is the final state of cycle 1. Firstly, we apply the control 
field for time period r 2 . Using (|24[) for m = 2 and assuming d 2 3> K.u; 2 is satisfied, we 
obtain the state 

\fo) =t/ (2) (r 2 )K) 

= cos(fi 2 T 2 )a 1 1) |1) +cos{n 2 T 2 )e- lE2T2 a i 2 ) |2) 

-isin(r! 2 r 2 )a( 1) |3). (31) 

We turn off the external field, then after free system evolution for time period r 2 , the 
final state is 

W 2 ) = e- iH ° T > |^ 2 > = 4 2) |1) + a 2 2) |2) + 4 2) |3) , (32) 

where 

af ] = cos(0 2 r 2 )e~ lBir 2a( 1) 
4 2) = cos(fi 2 r 2 )e- l£;2(T2+ ^ ) 4 1) 

4 2) = -isin{n 2 T 2 )e- lE ^a [ l ) (33) 
5.5. From (m — l)-f/i cyc/e io m-i/i cyc/e 

To figure out the explicit expression of Target state when all the control processes end. 
The recursion relation between coefficients of the (m — l)-th cycle and m-th cycle is 
need. Suppose after (m — l)-th cycle, we write down the state as 



,(m-l) 

fc=l 

Interacting with control field for time period r m and the restrictions al m 3> huj mi which 
lead to huj m /Q m — > 0, d m /fl m — > 1, we have the state 

\4>m) = COs(n m T m )a'f l ~ 1) |1) 

m 

+ ^cos(r» m r m )e-^ r ™4 m - 1) |fc) 

fc=2 

— i sin(f2 m T m )4 m 1 |m + 1) . (35) 
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After free evolution for time period r' m , the final state of the m-th cycle is 

m+l 

|^)=e-^|^) = ^4 m) |fe)- (36) 

where the coefficients are 

<4 m) = cos(r! m r m ) e - i; ^< ia ( m - 1) , (37) 

4 m) = cos(fi m r m )e- l£;fc ( T '"+ r ".)4 m_1) , 2<k<m, (38) 

a|T+i = - l sin(fi m r m )e-^+ 1 <"al m - 1) , (39) 
hereafter. Eqs. (l37H3"9"| establishes the relationship between the probability amplitudes 
of the (m — f )-th cycle and the m-th cycle. 

3.4- Target state 

The target state of 2 and 3 dimensional system has be given out in subsection B and 
C. The case N > 4 will be investigated in rest of this section. From |2"91 and |3"T1 we can 
find that 



a[' n ^ = Y\ cos(fliTi) exp 



(40) 



As o&_x J us ^ depends on a^ n , it can be written down as 



m — 1 



x exp 



i=l 
m — 1 



(41) 



For the coefficient using (|33p . we obtain 



Am) 



i sin(f2iTi) JJ^ cos(f^Ti 



i=2 



x exp < -1E2 



(42) 



To derive coefficient a[™ 1 ', 3 < fc < m, we use 



,(fe-i) 



isin(ft fe _iT fc _i) 



fc-2 



JT^ cos(r2iTi) exp 



fc-2 



(43) 



which is obtained from (141 1) by replacing m by fc — 1. Then we can recursively have 



cos(r2jTi) exp 



i—k 



-iE k 2J T, + t; 
= -isin($l fc _ir A; _i) cos(f^T 2 ) 

fc — 2 / m m N 



x exp 



\i=k i=k—l 



(44) 
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From expression (j33|, Notice that it does not recover the case k = 2, only to be valid 
for k = TO + 1. So far, all the probability amplitudes after m-th cycle are given by 1401 
[Handgl 

Since the system we considered here is N dimension, we implement N — 1 cycles 
to drive the system to arbitrary target states. Thus letting m = N — 1, we obtain the 
probability amplitude of the target state 



N-l 



(N-l) 



(N-l) 



(N-l) 



\ J cos(fliTi) exp 



N-l 



iEi 1 

i=l 

zsin(f2iTi) Y\ cos(fiiTj) exp < — iE 2 

i=2 { 

N-l 

isin(S! fc _i-r fe _i) cos(f2jTj) 



N-l 



N-l 
i=2 



i = l 



x exp 



/c-2 



'AT-1 



JV-1 



(3 < k < N). 



(45) 



5.5. Control parameters 



For a control problem, the target state, or in other words, the amplitude a^ _1 of the 
target state, is given. What we need to do is to determine the control parameters 
\ji,T[ I i = 1, 2, N — 1} from the probability amplitude of the target state. For 
convenience, we write the target state as 

N 

\^) = Y,lnC n \n) (46) 



71=1 



where C n s are the real part of the amplitude 

N-l 

Ci = JJ cos(rjiTi), 



N-l 



C 2 = sin(^iTi) Y[ cos(fiiTi), 

2V-1 

C„ = sin(ri fc -iT fc _i) ] I cos(f2iTi),3 < n < N, 



(47) 
(48) 
(49) 



i = l 



and 7„s are phases 

71 = exp 



N-l 



1E1 £ < 



i=l 

rjv-i 



72 = exp <; -iEi 



7„ = exp 



. i=2 

fc-2 



'AT-1 



jv-i \ 

E-i-? 



i—k 



i=k-l 
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(3 < n < N). (50) 

For a given target state, equivalently, C„ and 7„ are given, we can calculate control 
parameters {r„, r' n \n — 1, 2, N — 1}. From (|47 |) . (l48|) and Ci, C*2, we can determine 
T\ . Then form (l49l) with n = 3 , (|47|) and Ci , C3 , we can obtain t 2 . Repeating this 
process, we can find out all parameters r„, n = 1,2, N— 1, determining the coupling 
time between the system and field. 

All t[ can be obtained from (|50[) . From 71 and 72, we can obtain 

JV-l JV-l 

ExJ^rl, E 2 J2ri + E 2 ri, (51) 

i=l i=2 

from which we get the value of r{ and J^—j 1 t[. From 74, we find 

2 /JV-l N-l \ 

i=l \ i=4 i=3 / 

from which as well as r{ and ~^2,f = ^ t[, we can obtain t' 2 . Repeating this process, we 
can obtain all t[ i.e. all the time for free evolution. 

4. System II 

4-1- Time evolution operator 

For system II, the interaction operator is given in E.q.fJS]). This operator is same as 
that for system I except the state |1) is replaced by |m). So we can follow exactly the 
same procedure as in last section, namely, to simplify the time evolution operator 
In the m ~ th cycle, We define 

hu) m — E m+ i — E m (53) 

and reset the energy base point \{E m+ i + E m ) = 0. Using these two conditions in 
each cycle, we get the simplified free Hamilton: 

N 

H Q = ^E n \n)(n\ 

n 

= E m \m)(m\ + E m+ i \m + l)(m + 1| 

N 



^ m \ I I 1 ^ m I 1 1 \ / 111 

— - — \m){m\ H — \m + l){m, + 1 

N 



n = l 



TlUln 



N 



E E n \n)(n\ (54) 



n = 2 



and the interaction Hamilton is 

= d m (\m + l)(m\ + |m)(m + l|) = d m a^ n \ (55) 
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Thus the total Hamilton takes the form: 

JV 

-a r <+ £ E n \n)(n\+d m *W 



n = l 



N 



H '{m) + J- E n \n)(r 



n=l 



where 



(56) 



(57) 



With respect to the fact that the control field is constant, i.e. H dosen't contain 
t, the time evolution operator can be obtained as 

/ 

H '( m )+ E n \n)(n\\t 



U (m \t) 



exp 



V 



n = l 
n^n,, m + l 



cos{n m t) - -f - sin(fi m t)ffV) 



x exp 



JV 



-i £ E n \n)(n\t 



m > 2 



using the fact 



n 2n T 

Til ' 



. . 2n+l ,, . 



(58) 



(59) 



^.2. Determine amplitude a m 

From (|58l) for the case m = 1, [/^(i) is written as 



cos(Oi0 - sin(^it)^ (1) 



x exp 



JV 



:J2E n \n){n\t 



n—3 



(60) 



For this model, cycle 1 is exactly the same as the system I. So after the cycle 1, the 
system is driven to the state 

|Vi)=4 1) |l)+4 1) |2), (61) 

where 

<4 X) = cos(0 1 r 1 ) e - 4£lT S 

4 1} = -ism^n)^ 27 ^. (62) 

Different from system I, in cycle 2, the control field f(t) = di causes transition 
between |2) and |3). We can find the state after the cycle 2 in Schrodinger picture as 

IV^)=«i 2) |l)+4 2) |2)+4 2) |3), (63) 
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a^' — cos(f2 2 T2)e 



4 2) = cos(^ 2 r 2 )e- 4£;2 ^4 1) , 



= — i sin(f22T2)e 



-iE 3 r'(l) 
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(64) 
(65) 
(66) 



To obtain the target state, we first deduce the recursion relations between the 
(m — l)-th cycle and the m-th cycle. To this end, we suppose that, after m — 1 cycles, 
the system is on the state 



|C-i) = E4 m " 1) i^)- 



(67) 



k=l 



Then after interactions with the control field for time period r m , and free evolution 
for time period r' , we find the final state after cycle m as 



m+1 



with (m > 2) 



wj = E a i m} i fc ) 
fc=i 



ai m) = cos(fi m r m )e-^(^+^)ai" l - 1) , 1 < k < m - 1, 
a™ = cos(fi m r m )e-^ r ™a^- 1 ), 



(m) 



-ism(p, m T m )e 



-iE m+1 T' m {rn-l) 



(68) 



(69) 
(70) 
(71) 



From those recursion relations, and initial conditions (|62p we can get all the explicit 
expressions of . It is easy to see that 



Am) 



Am) 



Y\ cos(QiTi) exp < -iEi E( Tl + r-) 

i=l I L 

— i sin(iliTi) Y[ cos(fijTi) 



(72) 



x exp < -1E2 



J2(n + r[) + r' 2 + t[ 



i=3 



(73) 



= -»sin(r! m r m )e- lS ™+ 1 ^ a ^- 1 ) 



(— i) m 1 J^sin(f^Ti) x exp 

i=2 

m 

(-i) m sin(OjTj) exp 



-i E E i+i T i 

\i=2 

' m 

E^+i^ 



(74) 



Then using [6"9l and a^T , which is obtained from (|74p by replacing m by 771 — 1, we 
have 



a k J! cos (^i T «) ex P 



i=k+l 



iE k E (n + r[) 



i=k+l 



iE k r!Jk-l) 
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fc-i 



i=k+l 

/fe-1 



TiJ exp 



i=l 



vi=l 



fc-1 



(-i) 1 Y[ cos(fijTj) sin(OiTj' 



x exp < — iEk 



£( 

i=fe+l 



fc-i 



-IT, 



(75) 



One can check that (I75|) includes the case k = 2 and k = m &s special cases. 
Therefore, after N — 1 cycles, we arrive at the target state 



JV-l 



a[ N = Y\. C0S (^i T i) ex P { —iEi 



i=i 



N-l 



7V-1 



{-i) n 1 cos(^iTi) sin(fW 

i— n i— 1 

N-l 



x exp < — iE; 



£ 

,i— n+l 

(2 < m < N- 1), 

N-l 



n-l 



(AT-1) _ f _nN-l 
N 



(~i) N 1 Y\_ sin(fiiTi) exp 



'N-l 



(76) 



(77) 



4-3. Control parameters 

To determine control parameters Ti, t[, 1 < i < N — 1, we write the target state as 



AT 



JV 



IV') = E a « i n ) = E ^ nC n i n ) 



(78) 



in which 



N-l 



ra — 1 



C?l = TT COs(f2jTi 
i=l 
N-l 

C m = J| cosfe) Y[ wa(fliTi), (2 < m < N - 1), 

i— m i— 1 

N-l 

Cat = sin(^Tj), 



(79) 



i=i 



and phase j n 



rN-i 



71 = exp < — iE\ 



E^+ r «)+^ 



. i=1 
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/ -\n— 1 

(— i) exp < 



1 -\N-1 

In = (-») exp 



N-l 



m— 1 



i=i 



/JV-l 



. i=l 



(80) 



For a given target state, namely, C n and 7„ are given, we can determine the 
control parameters {r„, r^|n = l,2,...,iV — 1}. From Ci, C 2 we can determine ti , and 
then from C 2 , C 3 we obtained r 2 . Recursively we can obtain all r„. 

For r' n , from 71 and 72, we obtain 

JV-l JV-l 

E iJ2 T z+ E ir[ , E 2 J2< + E 2 t[ , (81) 



respectively. As E\ 7^ E 2 , we obtain t{ and 'r t '. From 73,74, we can obtain 



JV-l 



JV-l 



(82) 



i=3 



from which we obtain t 2 as well as X^i^s 1 T i- Repeating this process, we can obtain 
all t' . 



5. Conclusion 



In this paper we proposed a protocol to drive two types of finite dimensional quantum 
system to an arbitrary given target states. The control parameters are time periods 
{t to , T' m \m = 1,2, - ■ ■ ,N — 1} which can be explicitly determined from the probability 
amplitudes of the given target states. Relationship between control parameters and 
amplitudes is trigonometric functions and can be solved explicitly. 

We have 2{N — 1) real control parameters. In the target state there are N 
complex or 2N real parameters. Taking into account the normalization condition of 
target state, one has 2(N — 1) real parameters, the same as the number of the real 
control parameters. From this fact we can conclude that we can drive the system to 
an arbitrary target state by choosing appropriate control parameters {r m ,r^}. 

As further works, we would like to consider the indirect control protocol of finite 
quantum system by generalizing the control scheme in this paper. We also would like 
to consider the control protocol in the presentence of environment. 
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